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We study the properties of 74 Kr, and particularly the Gamow Teller strength distribution, using a 
deformed selfconsistent HF+RPA method with Skyrme type interactions. Results are presented for 
two density-dependent effective two-body interactions, including the dependence on deformation of 
the HF energy that exhibits two minima at close energies and distant deformations, one prolate and 
one oblate. We study the role of deformation, residual interaction, pairing and RPA correlations on 
the Gamow Teller strength distribution. Results on moments of inertia and gyromagnetic factors, 
as well as on EO and Ml transitions are also presented. 
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I. INTRODUCTION 



The study of exotic nuclei, characterized by nonoptimal N/Z ratios, has attracted an increasing attention during 
the last decades from both experimental and theoretical points of view jjJB]. Nevertheless, it is still true that most 
of the knowledge we have at present about nuclear systems corresponds to the nuclei inside the valley of stability, 
where the balance between protons and neutrons make these systems stable against particle emission. There is still 
much to be learned about systems out of this region of stability. Many complementary experimental projects and new 
facilities involving radioactive nuclear beams are being carried out or planned at CERN, GANIL, Grenoble-Munich, 
(T) ' GSI, Louvain-la-Neuve, RIKEN... to extend this limited knowledge. 

Information on these nuclei is of prime importance to understand different problems in nuclear physics related to 
the way nuclei respond as we vary their masses and isospin parameters to extreme values of the N/Z ratio. One can 
also learn about specific components of the N-N force and the modification of this force in the nuclear medium. Since 
the isospin dependence of the N-N interaction is largely unknown, the structure of single particle states, collective 
modes and behaviour of global nuclear properties is very uncertain in nuclei with unusual N/Z ratios. The theoretical 
techniques and parameters are optimized to reproduce the properties of known atomic masses and it is not clear how 
accurate is the theoretical extrapolation to regions far from stability. The study of nuclei far from stability provides 
c*j , then new basic pieces of information on nuclear structure. 

Knowledge of the decay properties of radioactive nuclei is also essential for nuclear astrophysics ||, where one 
tries to understand the production of energy and the synthesis of elements in stars during stellar events. In that 
respect, reactions with radioactive partners become very important, specially during violent stellar events with extreme 
conditions of density and temperature because the average time between successive nuclear reactions is much shorter 
than the average decay time of the radioactive nucleus. Therefore, these radioactive nuclei do not have time to decay 
before participating in new nuclear reactions. In particular, the decay properties of proton rich nuclei are fundamental 
to understand the r-p process (rapid proton capture nucleosynthesis) that involves multiple proton capture reactions 
on proton rich nuclei. The microscopic structure of exotic nuclei near the drip lines is therefore a topic of great interest 
in nuclear structure as well as in astrophysics. 

The neutron deficient side of the valley of stability is particularly interesting because in this region a large fraction 
of the Gamow Teller (GT) strength is accessible in (3 + decay j|. Usually kinematics limits j3 decay to a narrow energy 
window that in general misses the main part of the giant GT resonance, but this is not the case in proton rich nuclei. 
Therefore, proton rich nuclei provide a good opportunity to test the GT strength measured by other less direct means 
such as (p, n) reactions, where the GT strength is extracted by extrapolating the empirical proportionality between 
the zero degree cross section and the f3 decay strength to low- lying GT states || ■ 

All of these features, together with the ongoing experimental effort to measure GT strengths in proton rich medium 
mass nuclei ||, make this region the focus of this paper. 

It is clear that, at some point, extensions of the standard approaches will be necessary to cover particular aspects 
that, although not important in stable nuclei, may play a role in describing the unstable ones. Issues such as the role of 
the continuum Q in weakly bound systems or the neutron-proton pairing ||] in N rs Z nuclei should be considered. It is 
nevertheless important to know whether the most powerful tools, designed to reproduce the properties of nuclei within 
the valley of stability, are still valid when one approaches the drip lines. One should know the limits of applicability 
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of the existing interactions and methods. In our view this has not been sufficiently explored yet, particularly in what 
concerns the GT strength distribution. 

With this view in mind we study in this paper the (3 + decay of Kr, which is a medium mass proton rich nucleus 
of recent experimental ^ and theoretical interest. This nucleus is also representative of a mass region that exhibits 
a wide variety of shapes including large quadrupole deformations (fa — 0.4) and shape coexistence. Our method 
consists of a selfconsistent formalism based on a deformed Hartree-Fock (HF) mean field obtained with a Skyrme 
interaction including pairing correlations in the BCS approximation. We add to this mean field a spin-isospin residual 
interaction with a coupling strength derived by integrating over the nuclear volume the Landau-Migdal force, obtained 
from the same Skyrme interaction. The residual force is therefore consistent with the mean field. The equations of 
motion are solved in the proton-neutron quasiparticle random phase approximation (QRPA) Most (3 decay 

calculations in the past used as input in one way or another experimental single particle levels, but the extrapolation 
of these empirical information to exotic nuclei is at least questionable and the need for selfconsistent calculations 
in unstable nuclei has been emphasized (see for instance |ll]]). There are also selfconsistent mean field studies [[zj, 
as well as investigations of the dependence on deformation of the GT strength distributions in Tamm Dancoff 
approximation (TDA), but to our knowledge this is the first time that the selfconsistent deformed QRPA formalism 
described here is used to calculate GT strengths. A similar formalism for spin Ml strength distributions has been 



successfully tested |13|l4j] against experimental data in the rare earths and actinide nuclei. 

The paper is organized as follows: In Section 2 we introduce the theoretical framework used to describe the /? decay. 
We first introduce the selfconsistent method for the mean field and residual interactions and then we obtain the (3 
strengths in bare two quasiparticle (qp) approximation, in TDA, and in QRPA. In Section 3 we study the ground 
state properties and the low- lying structure of 74 Kr and compare our calculations to experiment. In Section 4 we 
discuss the results on Gamow Teller strength distributions and sum rules obtained in our calculations studying the 
sensitivity of these results to the various ingredients of our approach such as the effective two-body interaction, the 
pairing interaction, the nuclear shape, and the RPA correlations. We also show results on Ml transitions. Section 5 
contains half-life results in various approaches. Section 6 contains the main conclusions of this work. 



II. BRIEF DESCRIPTION OF THEORY AND DETAILS OF CALCULATIONS 



In this section we summarize the theory involved in the microscopic calculations presented in Sections 3 and 4. Our 
method is selfconsistent in the sense that, as described in Subsection A, with a given two-body density-dependent 
effective interaction we derive: a) the selfconsistent mean field generating the single particle energies, wave functions 
and occupations of the ground state, b) the particle- hole Landau-Migdal interaction averaged over the nuclear volume 
generating the QRPA modes. Thus, the same two-body interaction is used to derive the QRPA excitations consistently 
with the quasiparticle basis and ground state, as described in Subsection B. 



A. Mean field calculations and residual interactions 



It is well known that the density-dependent HF approximation gives a very good description of ground-state 
properties for both spherical and deformed nuclei |l5| and it is at present the most reliable mean field description. We 
consider in this paper two different Skyrme forces. On the one hand we use the Skyrme force Sk3 |ll| because it is the 
most extensively used Skyrme force and we consider it as a reference. We use Sk3 in its density dependent two-body 
version that has better spin-isospin properties than the three-body one [ [l3||l4j ]. On the other hand we use the SG2 
force |l7) of Van Giai and Sagawa. The two forces were designed to fit ground state properties of spherical nuclei 
and nuclear matter properties but, in addition, the force SG2 gives a good description of Gamow Teller excitations 
in spherical nuclei fT^ |. Recently |l3| ], we applied these two forces and a similar method to make a rather extensive 
study of isoscalar and isovector spin Ml excitations in deformed nuclei obtaining a good description of the available 
data, particularly with SG2. This gives us confidence on the predictive power of the method and on the reliability of 
the SG2 force when considering isovector spin properties of deformed nuclei. The parameters of these two interactions 
are given in Table 1, and the corresponding HF energy density functional for an even-even nucleus has the form 

£ (r) = ^ Pst ^2 | -*o/Vt' [1 - dss'Sw + xo (S ss > ~ 5 tt <)] 

st s't' 

[t s , v + Jv 2 p s 't' ) [1 + Sss>5 tt > + x 2 (S ss , + 5 tt >)) 
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+jxh {Ar a > v - 3V 2 p S 't') [1 - Sss'S tt < + x\ (6 SS < - 8 U >)] 
+-^hp a Ps't> [1 - $ss>Stt> + x 3 (S ss / - S tt ')} 

VoV-J 5 , t >(l + fe/))+£ c (r) (!) 
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with £c the Coulomb energy density 

-1 1/3 



-P P 0) 



(2) 



The spin-isospin (si) components of the nucleon, kinetic energy, and magnetization densities are 

^(rH^I^M)! 2 (3) 

i 

r s4 (r)=^^|V0»(r, S ,i)| 2 (4) 

i 

J st (r) = E «? # ( r > s '> *) H V x o-) & ( r > *. *) (5) 

with 



and similarly for r and J. 

In our calculations time-reversal and axial symmetry are assumed. The single-particle wave functions are expanded 
in terms of the eigenstates of an axially symmetric harmonic oscillator in cylindrical coordinates. The single-particle 
states \i) and their time reversed \i) are characterized by the eigenvalues fl of J z , parity TTi, and energy q 

l»)=E - 9 1 E CV„ 3AS |^ r n,AS) (8) 
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N n r ,n z ,A>0,'S 

with £li = A + £ > |, and 

|i)=E - 2 ^ C^ AS (-l)^ S |7Vn r n 2 -A-E> (9) 

iV n,.,n i ,A>0,S 

with n- t = — f7i = -A - £ < -|. For each iV the sum over n r ,n z ,A > is extended to the quantum numbers 
satisfying 2n r + n z + A = N. The sum over N goes from TV = to N = 10 in our calculations. 
The single-particle energies and wave functions are obtained from the HF equations 

wr e^ (10) 

We use the McMaster code that is based in the formalism developed in Ref. Ill] as described in Ref. (l9) . The method 
includes pairing between like nucleons in the BCS approximation with fixed gap parameters for protons, A p , and 
neutrons, A„. The number equation in the neutron sector reads 
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where vf are the occupation probabilities 



1 - 



E; 



uj = 1 



in terms of the quasiparticle energies 



Ei 



(12) 



(13) 



Eq. ( |ll| ) is solved iteratively at the end of each Hartree-Fock iteration to determine the Fermi level A„. Similar 
equations are used to determine the Fermi level and occupation probabilities for protons by changing TV into Z, A n 



into A p , and A„ into X p . 



The fixed gap parameters are determined phcnomcnologically from the odd-even mass differences through a sym- 
metric five term formula involving the experimental binding energies pCfl: 
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A n = -[B(N-2,Z)- AB (N - 1, Z) + 6B (N, Z) 

o 

-AB (N + 1, Z) + B (N + 2, Z)\ 



(14) 



A similar expression is found for the proton gap A p by changing N by Z and vice versa. For Kr we obtain 
A„ = A p = 1.5 MeV. 

The HF theory gives a single solution which is the Slater determinant of lowest energy. To allow for shape coexistence 
one has to extend the theory to a constrained HF theory |2lJ] . Minimization of the HF energy under the constraint of 
holding the nuclear deformation fixed is carried out over a range of deformations. When more than one local minimum 
occurs for the total energy as a function of deformation, shape coexistence results. The energy surfaces as a function 
of deformation that we will discuss in Section 3, are obtained by this procedure including a quadratic quadrupolc 
constraint [f2l|| . 

Following Bertsch and Tsai p2j the particle-hole interaction consistent with the HF mean field can be obtained as 



V. 



ph. 



16 ^ L 



-1) (T 1 • (T 2 



I + C-I^'t!-^ 



8 2 £ 



Sp st (ri) 8ps't> fo) 



(15) 



This gives a local interaction that can be put in the Landau-Migdal form |g3| . For the study of j3 decay the relevant 
residual interactions are the isospin contact forces generating the allowed Fermi transitions (AL = 0, AS = 0, AI* = 
0+) 

V F {12)= X F (tfta +*r*2) ( 16 ) 
and the spin-isospin contact forces generating the allowed Gamow Teller transitions (AL = 0, AS — 1, A/ 7r = 1 + ) 

V GT (12) = xgt o-i • <x 2 (t+tj + t^tt) (17) 



where we use the convention t + \p) 
spin-spin interaction Hss 



t \n) = \p) . The latter (Vgt) is the charge changing component of the 



Hss = [(1 + q) Si-s 2 + (l-q) s r s 2 ri • r 2 ] 



(18) 



considered in Refs. [|T^,0 for the study of spin Ml excitations. Not allowed transitions (AL > 0) produce strengths 
which are orders of magnitude smaller than the allowed ones (AL = 0) and will not be considered in this work. 

After functional differentiation we get from Eqs. (|]) and (]l5|-pT|), assuming symmetric uniform nuclear matter and 
averaging over the nuclear volume V 



XF 



AttR 3 



t (1 + 2x ) - -k% [t 2 (1 + 2a; 2 ) - h (1 + 2xi)] + -t 3 p a (1 + 2x 3 ) 
2 b 



(19) 



XGT 



to + ^4 (*i - 1 2 ) + ^hp a 



(20) 



4 



1/3 

where R is the nuclear radius and k>F the Fermi momentum kp — (3tt 2 p/2) . These coupling strengths are related 
to the familiar Landau-Migdal parameters F$ and G f (see for instance [O) by 

XF = m- ' XGT = vN- (21) 

where V = 4nR 3 /3 and No = (2TO*fcj?/?i 2 7r 2 ), with m* the effective mass. For compieteness the values of xf and 
Xgt are also given in Table 1. 

B. QRPA equations and j3 decay strengths 

Let us first consider Gamow Teller excitations. Using the notation = er^ (K — 0, ±1), we write Eq. ([i7]) in 
the more traditional form 

V GT = 2 XGT (-1)^ fePlK ( 22 ) 

K 

where in second quantization we can write 

/3+ = ^(n\a K \p) a+a+ = ^(n\a K \p) {u n v p a+a± 

np np 

+v n u p a n a p + u n u p a^a p + v n v p a>pa n } (23) 

and f3Z K = (—Tl) K (Pk) + ■ We use the standard notation a + (a) for particle creation (annihilation) and a + (a) for 
quasiparticle creation (annihilation) operators. We write the proton-neutron QRPA phonon operator for Gamow 
Teller excitations in even-even nuclei 

r ic = E 1 X "k< K< A ^\ (24) 

IK 

where A+ , Aj K are two-quasiparticle operators of the form 

= <<4 (25) 

Aj K = (-a^a+) + = a n a p (26) 

with fi n + fip = K, where n (p) stands for any neutron (proton) states \i) or |i) and h (p) is its time reverse. The 
QRPA equations for the even-even system are 

<0o \A 1K [H, F+ K ] I O ) = lok <0o \A 1K T+ K | 0q) (27) 

(0o | [H, T+ K ] A+ K | 4> Q ) = cj k <0o K K A+ K | fo) (28) 

with 

H = a tc*nE n + a+a p E p + V GT - (<f>a \H\ <f> ) (29) 

n p 

\(j) ) is the HF+BCS ground state satisfying a t \4>o) = for any t. From these equations the forward and backward 
amplitudes are obtained as 

X»« = 2X °1 {a lK Ml K + b 7K M"_") (30) 



Kk = 2 17 %kMI* + a 1K M"_«) (31) 

LOR T c 



5 



with £ lK = E n + E p the two-quasiparticle excitation energy with Ei given by Eq. (|l3|), A4± K are 

= E Kc*?£ + 6 7K >T/) ( 32 ) 



7-K 



- E + «7^ 7 T) (33) 

with 



7if 



and 



ij K = u n VpT, 7 ^ (34) 



(35) 



Ef = (n|(7K|p) (36) 



Owing to the transformation properties under time reversal and hermitian conjugation it is easy to show that: i) 
the K = — 1 and K = 1 modes are related to each other through time reversal and are degenerate; ii) the K = 
modes are time odd and satisfy = —X^° and = —Y"°. Forward and backward amplitudes are real. This 
allows to simplify the calculations to include only neutron-proton pairs with Q n > 0, i.e., the complete set of "/k 
neutron proton pairs can be reduced to the basic sets of pairs: 

• io — n.iPi with f2„ ; = f2 Pi > i 

u i_i = n.iPi with Q, ni — Q, Pi — 1 > ^ 

{riiPi with Q ni = Q Pi + 1 > | 
n-ipi with Q ni = Q Pi = i 



The matrix elements are 



Er<= J2 C^ AE+K C^ AE (2S)VTT^| (37) 



«\ = E ^Woi^oi ("^) (38) 



We therefore obtain the following dispersion relations for the if = and K = 1 modes: 



/l / o / ^ 2 C2 °*0 "T" I / , u io u io o C2 



for K = 0, and 



la ^o 



1 \ " V i " l ») c , I \ ■* 

o / / ^ 2 c2 V ' I / > ""V"'? 2 c 



a? 



V 5£ 5£— V — ^ !£— (40) 

Z^ ,,,, ,,,, _p. /_/ ,,,, + p, M , _r, v ; 



t p (p=±ij \ p p / jp(p=±i) \ p 



for if = 1. 



G 



The forward and backward amplitudes are then obtained from Eqs. (|3C|) and (|3lj) with the following normalization 
conditions 

2£[w) 2 -(>r) 2 ] = i ( 4i ) 



Using the inverse transformation 



E (K) ~{ Y C) = l ( 42 ) 

i p (p=±l) 



where stands for the time reverse of T UJK1 and the RPA conditions T U]K |0) = 0, \u)k) — ^t> K |Q) > w ith |0) the 



QRPA correlated ground state, one gets the /3£ strengths 



|/3± | 0) = (44) 

with 

M7=2^(a i0 X^+6 i0 F,7) (45) 



£ +w) (46) 

V(P=±1) 

A4" K is obtained from A^* exchanging aj K and &i K (see Eq. (0)). 

To finish this section we note that the Fermi transitions, although not shown here, can be handled in a similar way 
to the K = GT transitions with the important difference that the phonon operator for the Fermi mode is time even. 
Therefore, if we denote by ujp the frequencies of the Fermi excitations one has that X? F = X% F , Y? F = Y^ F and 
we may consider only pairs of the type io with the matrix elements a,i , bi replaced by 

a i = ( n i\Pi) (47) 
Kb = ( n i\Pi) V ni U pi (48) 



With these replacements the dispersion relation in Eq. ( |39| ) can also be used to find the frequencies u>f of the Fermi 
modes and the Fermi strengths can be obtained similarly to the K = GT strengths 

{co F |t+| 0) = MT = 2^ {a' i0 Xt; + b> i0 Y» F ) (49) 

< WF \t- 1 o> = m u j = 2^ k*t + «*) ( 5 °) 



with 



K: = KM? + b' i0 M«_ F ) (51) 



*r = - (*4^r + <m-_ f ) (52) 



and 
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X io I 



I io I 



with 



In Section 4 we show also for comparison the Gamow Teller strength distributions obtained in the TDA approxi- 
mation, as well as the uncorrelated two-quasiparticle neutron proton GT strengths. 

The equations for TDA are easily obtained from the above RPA equations by taking everywhere the limit 
1/ (uk + £i K ) — > 0, that in particular kills the Y amplitudes (Y?* = 0) and simplifies the eigenvalue equations. 
The strengths are then given by 

(u J T K DA \p%\^)=Mf DA (53) 

, TDA _ , , TDA 

M" + ° =2 1 £a i0 X£ (54) 

_ » / , TDA\ 2 

) =1 

io 

and corresponding relations for K = 1. 

Finally, the uncorrelated two-quasiparticle neutron-proton Gamow Teller (Fermi) excitations are obtained in the 
limit xgt = (\f = 0). Obviously in this case the excitation energies are the bare two-quasiparticle energies 



TDA 

M_° 



£ = £ lK = E n + E p , and the /J* strengths are 



{^k V \Pk\ <f>a) = a 
(u 2 ™\(3 K \fa)=b, 



(56) 
(57) 



Similar relations hold for Fermi transitions with K = and a 7o , 6 7o replaced by a' lo , b' lo (|47]) and ( |48| ) . 

In the uncorrelated case it is straightforward to show that the Ikeda sum rule is fulfilled. For Fermi transitions one 
has 



S- F ~S$ = I* I ^>o>| 2 - £ \(ur \t+\ 0o)| 2 = £ [Kf - «o) 



2 



E [K 



2 E K n< b* 



1 2 / 2 



«0 



For Gamow Teller transitions one has 

S GT — Sq t = £ (^GT,K ~ ^GT.Kj = ^ 
K=0,±1 

where we have used Eqs. (|56|) and (p7|) to write 



E<-E< 



A - 



3 (AT - Z) 



S'gT.K ~ S GT,K = E [K WA ' l^ 3 I ^o)) 2 ~ \( UK l^ + k°)| 2 ] = E [(V) 2 - ( a 7/c) 2 



(58) 



(59) 



(60) 



7 fi- 



It is a simple matter to show that the Ikeda sum rules hold also in the TDA and RPA approximations provided all 
the uj eigenvalues contained in the basis space are included in the sum so that the orthonormalization conditions 



El V^K Y UK V^K\rUK 
\ X 1K A 7'„ - Y 1K V, 



5^ Ktl ' in RPA 



and 



- , TDA TDA 

E x ™ x 7 = s ™,y k in TDA 



are satisfied. In practice, the strength functions are calculated up to an energy cut uj < E cut such that Ikeda's sum 
rule is satisfied up to a few per thousand. 
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III. GROUND STATE PROPERTIES AND LOW LYING STRUCTURE 



The constrained HF method allows one to get a solution for each value of the mass quadrupole Q-i- In Fig. 1 
we show the HF energy as a function of deformation for the two interactions SG2 and Sk3. The best HF solution 
at each Qi value can be obtained by varying the size and deformation parameters |D| of the deformed harmonic 
oscillator basis. As it is seen in Fig. 1 there are two minima, one in the prolate sector and one in the oblate sector, 
rather close in energy. Thus, as expected, shape coexistence takes place in this nucleus. This is an interesting feature 
that was previously predicted (see for instance Refs. |l]j24|] and references therein) and that has been experimentally 
confirmed |J. It should be noted that Fig. 1 leaves open the question as to whether the ground state is oblate or 
prolate. Sk3 favors prolate shape for the ground state and oblate for the shape isomer, in accordance with previous 
predictions J2^]. On the contrary SG2 force favors an oblate ground state and a prolate shape isomer. The results 
shown in Fig. 1 exemplify the situation met when several other Skyrme type interactions are used. In particular, the 
deformation energy curve obtained from calculations performed with the Skyrme force Ska [p6| shows a profile quite 
similar to that of SG2 with an absolute oblate minimum with Q 2 — —200 fm 2 and a prolate isomer at about 1 MeV 
with Q2 — 700 fm 2 . On the other hand the force Skm* [^7j produces a profile analogous to Sk3 with an absolute 
prolate minimum at about Q 2 = 700 fm 2 , and an oblate isomer at Q2 — —400 fm 2 . It is also worth mentioning that 
while the prolate solution appears at about the same deformation for all the forces considered, the oblate solution is 
spread from Q2 = —200 fm 2 up to Q2 = —500 fm 2 depending on the force considered. 

The experimental /3 2 values deduced from measured B{E2) strengths |2^] vary from |/3 2 | = 0.24 to | 1 = 0.40. 
These are consistent with our calculated parameter of deformation /3 2 defined in terms of the quadrupole moment Q2 
and the r.m.s. radius (r 2 ) 

The results obtained for /3 2 can be seen in Table 2. 

Although in Rcf. || it is stated that the ground state is prolate and the shape isomer is oblate, we find that exper- 
imental evidence supporting this assertion is still lacking and that it would be interesting to investigate further this 
question from the experimental side. Of course experimental evidence would require measuring the static quadrupole 
moments of the lowest 2 + states, however some hints can be obtained from a comparison to data of the theoretical 
low lying spectra and E0 transition rates obtained with these two interactions. To this end we compare in Fig. 2 
the experimental |2^] level spectra up to 12 MeV with those obtained for each interaction neglecting shape mixing. 
Clearly, the low lying spectra in Fig. 2 obtained with the prolate Sk3 compare better with experiment. We follow the 
approximate variation after projection method ]3C[ | and write the energy of the K — + band head (J = + ) as 

/J2\ 

E 0+ =E HF -^—L (62) 

L.L CT 

and the energies of the higher (/ = 2, 4, 6, ...) members of the rotational band are calculated to lowest order in angular 
momentum |31| 

E I+ - E 0+ = (63) 

where T cr is the cranking moment of inertia, 

(64, 

with \atPt) = c-ti t a ^3 t I0o) any two quasiparticle state (t = proton, neutron). The values of the Fermi level for neutrons 
A n and protons X p , the values of (J 2 ) and I cr , together with the charge r.m.s. radii, deformations /3 2 , HF energies 
(Ehf) , and + energies (E +, from Eq. ((6^)), are given in Table 2 for the two minima, oblate and prolate, with the 
two interactions, SG2 and Sk3. The experimental (2^] binding energy is i?( 74 Kr) = —631.281 MeV. 

In any case, since the two minima are close in energy and the barriers are not so high, substantial shape mixing 
can be expected. In Fig. 3 we show that the experimental |9| p 2 (E0) value is compatible with both Sk3 and SG2 
predictions when shape mixing is taken into account. Considering that the |0^~) and [Oj) states mix to form the 
observed + ground state |0+) and first + excited state |0^ 7 ) 
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|0+) = VA|0+) + \/l — A|0+) 



(65) 



\o+j) = VT=\\of) - V\\o£) 

The amplitude of the monopole transition is given by pH 

1 



p (EO, 



3R 2 



01 



EO 

I 

- (2A- i) (o| #o Of 



VA(1-A)((0+ 



#o 



(66) 



of>- 



0+ 



#o 



(67) 



with £0 = ejf|. The crossed term in Eq. (|6j) is strictly zero with either of the forces Sk3 or SG2 in the HF limit 
(y 2 — 0, l) and can be neglected because of the large (3— spacing in between the Of and Of minima in Fig. 1. One 
therefore gets for the p 2 (EO) strength 



p 2 (EO) = A (1 — A) 



Z 
R 2 



1 2 



(68) 



with n and r-i the charge r.m.s. radii of the Of and Of states. In principle, R should be the r.m.s. radius of the ground 
state but since to our knowledge this has not been measured, we use R = (r\ + T%)l1 in Eq. d6§|). An approximate 
expression for p 2 (EO) in terms of /?2(i) an d /?2(2) na s been used in Ref. 



p 2 (£0)=A(1-A) 



2(1) 



$2(2) 



5^5 
2lV3 



2(1) 



^2(2; 



(69) 



In Fig. 3 we plot p (EO) as a function of A, obtained from both Eq. ( p8[ ) and Eq. (|69j), using the r.m.s. 
and deformation values obtained with either Sk3 or SG2 forces (see Table 2). One sees that using Eq. (^38j) for 
Sk3 the experimental EO value calls for a large mixing (A (A — 1) ~ 0.25) while for SG2 the mixing must be small 
(A (A — 1) ~ 0.05). This is consistent with the fact that with SG2 the two minima differ more in energy and the 
barrier is higher than with Sk3 (see Fig. 1), hence a much smaller mixing probability is to be expected. 

We would like to point out that also the Ml strength between the lowest 2 + states can be obtained in an analogous 
way from the gyromagnetic ratios gn of the Of and Of bands given in Table 2. The B(M1) strength between the 
lowest 2+ states would be [BlJ 



B (Ml, 2 
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47T 



A(l-A)[ 5 



R,l 



l 2 2 
9R,2\ Mat 



For Sk3 with A = 0.5 and g RA = 0.487, g R „2 = 0.459 (see Table 2), we get B (Ml, 2+ 
SG2 with A = 0.05 and g R>l = 0.498, g Ra = 0.463 (see Table 2), we get B (Ml, 2+ -> 2+ 



(70) 

2+) = 0.010 p 2 N and for 
= 0.003 p 2 N . 



IV. GAMOW TELLER STRENGTHS AND SUM RULES 



The j3 + strength distributions calculated in the selfconsistent HF+RPA scheme described in Section 2 are shown in 
Fig. 4 as a function of the excitation energy of the daughter nucleus. We present the results for both of the interactions 
considered, SG2 and Sk3, corresponding to /3 + decay from each of the two minima, oblate and prolate. For comparison 
we also present the results that would correspond to f3 + decay from the spherical shape. A strong shape dependence 
is observed in Fig. 4 that is roughly the same independent on whether the SG2 or the Sk3 interaction is used. The 
Qec values are indicated by vertical lines. The Qp value for (3 + decay has been calculated for each of the prolate, 
oblate, and spherical HF+BCS solutions as 

Qp+ = M p - M n - m e + X p - A„ - E p - E n (71) 

with A p (A n ) the proton (neutron) Fermi energy given in Table 2 and E p (E n ) the lowest proton (neutron) quasiparticle 
energy. The associated Q— value for electron capture is 
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Qec = Q/3+ + 2m e 



(72) 



The values of Qec obtained and shown in Fig. 4 are the following: Qec — 4.848 MeV for the SG2 oblate solution, 
Qec — 4.628 MeV for the SG2 prolate solution, Qec = 4.892 MeV for the Sk3 oblate solution, and Qec = 4.308 
MeV for the Sk3 prolate solution. It is important to stress that the strength summed up to energies below the Qec 
window is also strongly dependent on deformation (see also Fig. 8) and this should also be experimentally observable. 

Fig. 4 summarizes the main results of this section. One should keep in mind that a quenching of the §a factor 
(<7A,e// = 0.7 g a) is expected on the basis of the observed quenching in charge exchange reactions and spin Ml 
transitions in stable nuclei, where g s ,eff is known to be approximately 0.7g s j ree . 

In what follows we disentangle the meaning of these results in the light of simpler approximations. In particular, we 
will study the sensitivity of the Gamow Teller strength distribution to RPA correlations and to pairing correlations, 
as well as to the nuclear shape. 

Since the distinct profiles of the strength distribution in the three nuclear shapes considered are not seriously 
modified from one interaction to another we restrict the discussion in what follows to results obtained with the SG2 
interaction. 

In the laboratory frame the transition probability for /3 + decay from the + to a 1+ state is given by 

^T(0 + -li) = f||(li||/3 + ||0 + )| 2 (73) 

and the reduced matrix element (1+ + ) in the laboratory frame is related to the intrinsic matrix elements 
(lok \P~k \ 0) discussed in the previous section by 

|(C||/3+||0+)| 2 =2|< Wl |/3+|0)| 2 , w = «i 

= |(c;o|/3 + |0)| 2 , lu = uj (74) 

where we have used the Bohr-Mottelson factorization approximation |n| and have neglected possible higher order 
corrections due to angular momentum projection |3(J. Therefore, the calculated strength functions can be written as 

B+ T H = |||^|(c,o|/3 o + |0)| 2 ( 5(c,-c, o ) + 2^|(c, 1 |/?+|0)| 2 ( 5( W - Wl )| (75) 

L WO wi ) 

consisting of spikes at the various ojk {K = 0,1) solutions. 

To facilitate the comparison among the various approximations discussed here we have folded the calculated GT 
strengths with Y = 1 MeV width Gaussians converting the discrete spectrum into a continuous profile. This can also 
be interpreted as an approximate way to incorporate the further fragmentation and smoothing effects due to coupling 
to other excitation modes not taken into account in RPA. The effect of this folding procedure is illustrated in Fig. 5 
where we compare the strength in the discrete spectrum to the convoluted strength distribution for the particular case 
of the prolate solution with the SG2 interaction. In the lower panels the K = and K — 1 strengths are represented 
by vertical dotted and solid lines, respectively. The total strength in the upper parts of Fig. 5 is represented by a 
solid line. The right panels contain the results of the RPA calculations while the left panels contain the results of 
calculations for f3 + decay to uncorrelated neutron proton two quasiparticle (2qp) excitations. The large difference 
observed between the results of RPA and uncorrelated 2qp neutron-proton excitations in the discrete spectra in Fig. 
5 is reliably reflected in the folded distribution. Hence, from here on we restrict our presentation to folded strength 
distributions. 



A. Comparison of RPA, TDA, and bare 2qp strengths 



The importance of the role played by the residual interaction is illustrated in Fig. 6 where we compare results for GT 
strength distributions obtained in RPA (solid line), in TDA (dashed line), and in the uncorrelated two-quasiparticle 
limit (dotted line). Clearly the repulsive spin-isospin residual interaction (Vqt) moves the strength to higher energy. 
This effect is already present when one goes from the bare 2qp limit to the Tamm Dancoff approximation. The 
deformation dependence of the strength agrees with that found by Hamamoto and collaborators p2 j in the context of 
TDA. The main effect of the RPA correlations, that can be seen comparing TDA to RPA, is to provide a reduction 
of the TDA strength maintaining the position of the peaks practically unchanged. A strong dependence of the 
GT strength on the shape of the f3 + parent is observed in the three approximations, being more dramatic in the 
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uncorrelated situation. This reflects the fact that the dependence on the different internal shell structure involved in 
the different shapes is stronger in the absence of a residual interaction. The latter tends to redistribute the strength 
and tends to produce a smoother strength distribution. This is further illustrated in Fig. 7 where we show the RPA 
results for oblate, spherical, and prolate shapes corresponding to different coupling strengths of the Vqt interaction. 
The results corresponding to xgt — 0.48 in Table 1 are compared to those obtained when the strength xgt of the 
interaction is reduced or increased by a factor of 2. 

The summed strengths up to different excitation energies are represented in Fig. 8 for the two minima (oblate and 
prolate) corresponding to the SG2 interaction. The results with Sk3 are similar. Also plotted in the figure are the 
linear energy- weighted sum strengths. Clearly the uncorrelated sum strength 

E( s x T+ U orr = f|EK^k* + l^>| 2 (76) 

x a/3 

is conserved by the TDA approximation while the RPA approximation conserves the uncorrelated LEWSR 

E (3A OT + )un C orr = ^ E & a + E ) |(o/9 | «t+ \ ^ (77) 
x a/3 

where a(3 represents any two quasiparticle neutron-proton excitation. Although this is to be expected from general 
theory, it is interesting to sec that already at low energies (E* < 10 MeV) in RPA the linear energy-weighted sum 
tends to converge to the bare 2qp result, while the non energy- weighted sum is clearly lower than the bare 2qp result 
in order to compensate for the push to higher energies of the strength produced by the Vqt interaction. The opposite 
happens with TDA. The lack of RPA correlations makes the TDA linear energy-weighted sum to start diverging 
already at 10 MeV. Also shown in Fig. 8 and indicated by arrows, are the calculated Qec values for the oblate and 
prolate shapes where one sees that the predicted summed (3 + strength below Qec is larger in the prolate case. The 
total linear energy-weighted and non energy-weighted sums up to an excitation energy of 30 MeV are given in Table 
3. The results correspond to f3 + decay from the two minima (prolate and oblate) with the SG2 interaction. Also 
given for comparison are the results for (3 + decay from the prolate minimum with Sk3 interaction. Here one can also 
see (see Fig. 4) that the results depend more on the shape of the parent nucleus than on the effective interaction. For 
completeness in Table 3 we give also the (3~ summed strength and in the column labelled Ikeda, we give the difference 
between (3~ and (3 + summed strengths to show that the Ikeda sum rule is satisfied in all cases with an accuracy of 
0.3%. The summed strengths for K = and K = 1 are given separately to check that Ikeda's sum rule is satisfied 
independently for each K value. 

B. Analysis of pairing and deformation effects 

As already stated our mean field calculations only include neutron-neutron and proton-proton pairing correlations. 
Taking into account neutron-proton pairing would increase the diffuseness of the Fermi surface (see Goodman pi), 
which is governed by the gap parameters. It is therefore interesting to study the sensitivity of the GT strength to 
this diffuseness. To this end we compare in Fig. 9 the RPA results obtained with the SG2 force for different values 
of the gap parameters A„ = A p = 1.5 MeV, A„ = A p = 1.0 MeV, and A„ = A p = 0.5 MeV. In the spherical case 
the role of pairing is different in the low (E* ~ 1 MeV) and high (E* ~ 6 MeV) peaks. The first peak decreases with 
increasing pairing while the second peak increases. In the deformed cases (oblate and prolate) the role of pairing is 
somewhat less pronounced, but as a general rule, we see that with increasing pairing the strength at high energies 
(beyond E* ~ 4 MeV) increases. 

To understand better the role of pairing and deformation we show in Fig. 10 the results of bare two quasiparticle 
excitations with and without pairing for the various shapes. From left to right we can see the effect of going from 
prolate to spherical and to oblate. From down to upper panels we can see the effect of going from no pairing (A = 0) to 
pairing (A = 1.5 MeV). In these figures the solid lines (both in the discrete and continuous distributions) correspond 
to K = 1 excitations and the dotted lines to K = excitations. When there is no pairing and no deformation 
(Fig. lOe) there is only one neutron-proton particle-hole channel open (71753/2 — * wpxW- When we include pairing 
maintaining the spherical shape (Fig. 10b), the strength of this channel decreases and new channels are open, in 
particular the Trfv/2 —> vjhli that is forbidden in the case A = 0,/3 = (see also Fig. 11). Consequently a two peak 
structure appears and it remains in the spherical RPA results (see for instance Fig. 4), even though it is modulated 
by the effect of the residual interaction and RPA correlations. Deformation alone also causes the opening of new 
channels and the reduction of the 71773/2 ^ vpi/2 peak (see the transition from panel (e) to panels (d) and (f) in Fig. 
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10), but in addition it causes fragmentation of the strength between single particle spherical orbitals. In particular, 
the K = and K = 1 modes are degenerate in the spherical case independently of whether there is pairing or not. 
This degeneracy is clearly broken in the deformed case: only K~- — > transitions with = K v = 1/2 contribute 
to both K — and K = 1 modes, while — > Kf^ transitions with K v — K„ ± 1 (K v , > 1/2) contribute only to 
K = 1 modes, and those with _K"„ = 7^ > 1/2 contribute only to K = modes. Clearly which are the new channels 
that are open with deformation, and what are their strengths, depends strongly on whether the nucleus is oblate or 
prolate and on the magnitude of the deformation. As seen in Fig. lOd (see also Fig. 11) in the prolate A = case 
the strongest new channel is the K = 1 (7/2~ — > 5/2~), that has a strong spherical component 7r/ 7 / 2 — > vfzii-, while 
the 7rp 3 / 2 — > transition is negligibly small, and other positive parity states have also more strength. On the 

contrary in the oblate case the strength is still due to the 71753/2 — > vpi/2 transition that appears strongly fragmented 
into several 3/2~ — > l/2~ and l/2~ — > l/2~ transitions. The effect of pairing in the deformed case is similar to that 
in the spherical case. It opens more channels and reduces the strength of the existing ones (compare panels (a) to 
(d) and (c) to (f) in Fig. 10), resulting in more complex and smoother strength distributions. It is also interesting to 
observe that in spite of the redistributions of strength caused by the residual interaction Vqt and RPA correlations, 
the distinctive prolate, spherical, and oblate profiles of the Bqt strengths in Fig. 4 are still reminiscent of those in 
Fig. 10. Hence, the above discussion helps to trace back the origin of the main features observed in Fig. 4. 

In order to understand better the mechanisms leading to such patterns in the profiles of the GT strength distribution 
it is illustrative to compare the HF single particle neutron and proton spectra around the Fermi levels in the spherical 
and deformed cases. This can be seen in Fig. 11, where the Fermi energies are indicated by the step dashed lines. The 
labels correspond to the £j values in the spherical case and to K n in the deformed cases. In the absence of pairing, the 
allowed GT + transitions connect only protons below X p with neutrons above A„ with AK = 0,±1 and same parity. 
The most important GT transitions are indicated by the solid arrows. They are (71753/2 — ► ^Pi/2) m the spherical 
case, (7r7/2~ — * v5/2~) in the prolate case, and (7iT/2~ — > i/l/2~) in the oblate one. The dashed arrow indicates the 
most important GT transition in the spherical case —* v /s/2), which is allowed once pairing correlations are 

included. From this figure one can also understand why the GT strength in the oblate case is much smaller than in 
the prolate one. As seen in Fig. 11 the reason for this behaviour is that in the oblate case the proton levels below the 
Fermi energy are negative parity states while those above the neutron Fermi level are mostly positive parity states, 
while in the prolate case one finds both, positive and negative parity states, below X p and above A„. This figure can 
also be used to predict behaviours of GT strengths when we add or subtract nucleons to 74 Kr. For instance, when 
we remove two neutrons, the state v3/2~ in the oblate case will appear above the Fermi level and therefore it can be 
fed by (3 + decay of the negative parity proton states below the proton Fermi level. Thus, one expects an important 
increase of the GT strength when going from the oblate 74 Kr to the oblate 72 Kr. On the contrary, based on the same 
picture, we do not expect much changes in the prolate case. 



C. Magnetic properties 

Gyromagnetic factors and estimates of the total Ml transition strength between the lowest 2 + states were given at 
the end of Section 3. In this Section we focus on + — > 1 + spin excitations, which are stronger. 

The spin Ml transitions are the AT Z = isospin counterparts of the AT Z = ±1 Gamow Teller transitions. The 
study carried out in jl3| for several isotope chains showed that the Ml strength depends on deformation, not only 
for orbital but also for spin excitations. Thus, on general grounds, one may argue that the deformation dependence 
of the GT strength discussed earlier should be comparable to that of the spin Ml strength considered in (l3|] . 

To see to what extent this argument holds we show in Fig. 12 the profiles of the spin Ml strength distributions for 
the three shapes oblate, spherical, and prolate in 74 Kr. The results correspond to selfconsistent HF+RPA calculations 
with the SG2 interaction. These calculations are described in |l3|,[l4| and closely follow the description of calculations 
in Section 3 except for the AT Z = character of the Ml operator that, in particular, demands orthogonality of the 
RPA states to the spurious rotational state |Q. For completeness we also show the results obtained in the absence 
of residual spin-spin interaction (Kg = in Eq. (|l8|)), i.e., the results for bare two quasiparticle spin Ml excitations. 

One can see in Fig. 12 that as in the case of GT excitations (compare with Fig. 6) the residual interaction 
produces a redistribution of the strength pushing it to higher energies. There is a striking similarity between the 
shape dependence in Fig. 12 and that in Fig. 6. This is so in both 2qp and RPA strength distributions. 

Although it is unlikely that Ml transitions may be observed in such highly unstable nuclei, the comparison of Figs. 
6 and 12 tells us that the main features of GT and spin Ml strength distributions are very similar. This suggests 
that we may learn about properties that are observable in highly unstable nuclei (like (3 decay) from properties that 
are observable in stable nuclei (like Mi's), and vice versa. 



13 



V. TOTAL HALF-LIFE 



Another quantity of interest in nucleosynthesis is the (3 + / EC half- lives of exotic proton rich nuclei. The total 
half-life depends on the strength distribution in a different way from both the energy-weighted and non-weighted 
sums considered above. It is therefore illustrative to see the predictions for the half-life as well as the sensitivity of 
this quantity to the various effects considered in this work. 

The partial half- life t\i 2 for a given transition connecting the ground state of a nucleus with an excited state in the 
daughter nucleus is given by 

/ (Z, u) t l/2 = — — - -5 (78) 

where we have neglected the Fermi strength. / (Z,u>) is the Fermi integral (|]34|]. We use D = 6200 s and include 
effective factors 



K 2 



2 



(9A/gv) eff = 0.7(9A/9v) free = 0.74 (79) 



The total half- life T X j 2 for allowed f3 decay from the ground state of the parent nucleus is given by summing over 
all the final states involved in the process 

r i7> = ^£/(^)l<tfll^ll 0+ >l a (80) 

UJ 

The Fermi integrals / (Z, u>) are taken from Ref. P4| . 

One should notice that the sum over final states in Eq. ( |S0| ) is extended over all the states in the daughter nucleus 
with excitation energies below the Q— value. Another difference with respect to the energy-weighted and non- weighted 
sums considered earlier is that in the sum of Eq. ( |80| ) the weighting factor is larger at low excitation energies, which 
correspond to large kinetic energies of the emitted positron. 

We can see in Table 4 the results obtained from bare 2qp, TDA, and RPA calculations for the total j3 + /EC half-life 
of 74 Kr. Results are shown for the two Skyrme forces Sk3 and SG2, as well as for the two shapes oblate and prolate. 
The experimental value is (li/ 2 ) EXP — 11-5 min. It is clear that the half-lives calculated with both Sk3 and SG2 
increase going from 2qp to RPA and are larger in the prolate case. On the other hand, the oblate SG2 produces larger 
Ti/2 values than the oblate Sk3, while the prolate SG2 gives smaller Ti/ 2 values than the corresponding prolate Sk3. 

Although not shown in the table we have also checked the sensitivity of this quantity to other effects considered 
in this work such as the strength of the residual interaction xgt and the pairing gaps. For example, from an RPA 
calculation with the force SG2 we find that changing the strength xgt in Table 1 to half or twice its value, Txj 2 
changes from 3.5 to 16.7 min in the oblate case and from 6.9 to 23.8 min in the prolate case, i.e., roughly speaking, 
within RPA, Tyi 2 increases proportionally to xgt- The results are also sensitive to the pairing gaps. Again, an RPA 
calculation with SG2 and pairing gaps A„ = A p — 1 MeV, gives T\i 2 = 1.8 min and Ti/ 2 = 4.9 min in the oblate and 
prolate cases, respectively, i.e., the half-life decreases with decreasing pairing. 

One should also keep in mind that values for the effective qa factors ranging from 0.7 to 0.8 of the bare values are 
considered in the literature. Since we have taken here the value 0.7, the half-lives in Table 4 would be reduced a 20% 
by using the value 0.8. 



VI. SUMMARY AND FINAL REMARKS 



We have investigated shape isomerism and j3 decay in 74 Kr on the basis of the sclfconsistent HF+RPA framework 
with Skyrme forces. This is a well founded method that has been successfully used to describe quite diverse properties 
of stable spherical and deformed nuclei over the periodic table. The merits of this approximation are well known fl3|| - 
p2[ . Once the effective interaction parameters are determined by fits to global properties in spherical nuclei over 
the periodic table and the gap parameters of the pairing force are specified, there are no free parameters left. Both 
the residual interaction and the mean field are consistently obtained from the same two-body interaction. This is 
particularly desirable for nuclei far from stability, where extrapolation of methods based on phcnomcnological mean 
fields are more doubtful. With the method used here one predicts ground state and low-lying properties, as well 
as (3 decay and spin Ml excitations without introducing any free parameter. We obtain a fair description of the 
experimentally known properties of Kr and, in particular, of the interesting shape isomerism. In this regard one 
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question is open as to whether the ground state is prolate and the first excited + state is oblate or vice versa, since 
different Skyrme interactions give different answers. All the Skyrme interactions that we tried (Sk3, SG2, Ska, and 
Skm*) give a similar stable prolate shape (fo — 0.4), but their predictions on the stable oblate shape are somewhat 
different ranging from (3 ~ —0.15 up to (3 ~ —0.30. These deformations are compatible with the experimental B(E2) 
PH and p 2 (E0) || values. Theoretical predictions for gyromagnetic ratios and Ml transitions are also given. 

The RPA Gamow Teller (3 + strength distributions depends strongly on the shape (prolate, spherical or oblate) 
of the f3 + parent 74 Kr, and remarkably, this dependence is quite noticeable in the experimentally accessible energy 
window. A close similarity is found between this deformation dependence and the deformation dependence of the 
spin Ml strength distribution. It is important to emphasize that these results do not depend much on which effective 
Skyrme interaction (Sk3 or SG2) is used. 

We have also studied the dependence of the results on various aspects of the theory, in particular on the residual 
interaction, RPA, and pairing correlations. Compared to the uncorrelated two quasiparticle response, RPA shifts 
the GT strength to higher energies and reduces the total strength. These effects are of course more pronounced 
if one increases the coupling strength of the residual interaction. While the shifting effect is already contained in 
the TDA description, the quenching effect is not. This is consistent with the fact that the non energy-weighted 
summed strength is conserved in TDA, while the linear energy weighted summed strength is conserved in RPA. BCS 
correlations reduce the strength of allowed particle-hole excitations and create new ones; the main effect of increasing 
the Fermi diffuseness is to smooth out the profile of the GT f3 + strength distribution increasing the strength at high 
energies. 

Predictions for the half-life are given together with a discussion of the sensitivity of this quantity to the various 
effects considered in this work. The experimental value, (T 1 / 2 ) £ , xp = 11.5 min, is in between the two RPA results 
for the prolate solutions of the SG2 and Sk3 forces. 

It will be interesting to test our results with experimentally accessible information on GT strengths. Since results 
on the proton rich odd- A Kr isotopes will be soon available , it will also be interesting to extend our calculations to 
include odd-A nuclei. Comparison to data would be desirable before considering further theoretical refinements. The 
latter would deal with issues such as extended RPA, or the inclusion of the continuum and neutron-proton pairing 
at the mean field level. It is however important to keep in mind that in the present method the Ikeda sum rule is 
conserved, while violations of this sum rule are found in some of these extensions. 
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Figure 1. Total energy as a function of the mass quadrupole moment obtained from a constraint HF+BCS 
calculation with the Skyrme interaction SG2 (solid) and Sk3 (dashed). 

Figure 2. Experimental spectrum for 74 Kr p9| compared to the theoretical rotational spectra obtained for the 
oblate and prolate solutions of the SG2 and Sk3 interactions. 

Figure 3. EO strength from Eq. ( pq ) for SG2 (solid line) and for Sk3 (dotted line). Also shown is the EO strength 
from Eq. (|69| ) for SG2 (dashed linejand for Sk3 (dash-dotted line). The strength is plotted as a function of the 
mixing parameter A (see text). The experimental value from is also indicated. 



Figure 4. Gamow Teller strength distribution in 74 Kr as a function of the excitation energy of the daughter 
nucleus. The results correspond to the forces SG2 (solid lines) and Sk3 (dashed lines) in RPA with the coupling- 
strengths given in Table 1. The results are for the prolate (upper part), spherical (middle part), and oblate (lower 
part) shapes. The vertical solid (dashed) lines indicate the SG2 (Sk3) Qec values. 

Figure 5. Energy distribution of the Gamow Teller strength for the prolate solution of the force SG2 in 74 Kr. The 
lower panels show the spectrum of K = (dashed vertical lines) and K = 1 (solid vertical lines) excitations. The 
upper panels show the corresponding folded strength distributions (see text) for K — (short dashed line), K = 1 
(long dashed line), and total (solid line) GT strength. The left (right) panels are from bare two quasiparticle (RPA) 
calculations. 



Figure 6. Comparison of RPA (solid line), TDA (dashed line), and bare two quasiparticle (dotted line) GT 
strength distributions in 74 Kr. The results correspond to the force SG2 for the prolate, spherical, and oblate solutions. 
Also shown with vertical lines arc the Qec values. 



Figure 7. RPA Gamow Teller strength distributions in Kr for different coupling strengths of the residual 
interaction xgt- The solid lines correspond to the value in Table 1. The short (long) dashed lines correspond to a 
half (double) value of the coupling strength. 

Figure 8. Accumulated GT strength (Y.B [g^/in]) and energy-weighted strength (Y.EB [MeV g^/An]) in 74 Kr 
as a function of the excitation energy of the daughter nucleus. The results correspond to the prolate (upper part) and 
oblate (lower part) solutions of SG2. Dotted lines are for bare 2qp calculations, dashed lines for TDA calculations, 
and solid lines for RPA calculations. The vertical lines indicate the Qec values. 



Figure 9. Pairing effect in the RPA Gamow Teller strength distribution in 74 Kr. The solid lines correspond to 
calculations with pairing gaps A„ = A p = 1.5 MeV, long (short) dashed lines are for A„ = A p = 1.0 MeV (A„ = A p = 
0.5 MeV). As in previous figures we can see the results for the prolate (upper part), spherical (middle part), and oblate 
(lower part) solutions. 

Figure 10. Comparison of pairing and deformation effects in the Gamow Teller strength distributions of 74 Kr. The 
results correspond to bare two quasiparticle calculations with the force SG2. From left to right we can see the prolate, 
spherical, and oblate cases. The upper panels correspond to results with pairing (A„ = A p = 1.5MeV) and the lower 
ones are without pairing (A = 0). Solid lines (spectra and folded distributions) correspond to K = 1 excitations, while 
dashed lines (spectra and folded distributions) correspond to if = excitations. The main configurations leading to 
the GT excitations are also shown. 
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Figure 11. Spherical, prolate, and oblate single particle energies for protons and neutrons obtained from HF 
calculations with the force SG2. The Fermi energies are indicated by the step dashed lines. The most important GT 
transitions in the spherical and deformed cases are indicated by the solid arrows while the dashed arrow indicates the 
most important GT transition in the spherical case allowed once pairing correlations are included. 

Figure 12. Energy distribution of the spin Ml strength with K n = 1 + in 74 Kr. The results correspond to the 
force SG2 and pairing gap parameters A„ = A p = 1.5 McV. Dashed lines are for 2qp calculations and solid lines for 
RPA calculations. The three panels are for the prolate (upper part), spherical (middle part), and oblate (lower part) 
solutions. 
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Table 1. Parameters of the Skyrme forces SG2 and Sk3: t [MeV fm 3 ] , h [MeV fm 5 ] , t 2 [MeV fm 5 ] , t 3 [MeV fm 6 ] , 
W [MeV fm 5 ] , xq, x%, x 2 , %3, and a. Also shown are the strengths of the separable isospin xf [MeV] and spin-isospin 
Xgt [MeV] residual interactions obtained from Eqs.((l^) and ( po|) , respectively. 

to h t 2 t-j, W xo xi x 2 x 3 1/q xf Xgt 

SG2 -2645.0 340.0 -41.9 15595.0 105.0 0.09 -0.0588 1.425 0.06044 6.0 0.69 0.48 
Sk3 -1128.75 395.0 -95.0 14000.0 120.0 0.45 0.0 0.0 1.0 1.0 0.88 0.46 



Table 2. Fermi energies [MeV] for neutrons (A„) and protons (A p ),values of (J 2 ), cranking moments of inertia I cr , 
gyromagnetic ratios gn, charge radii r c , deformation parameters /?2, total energies Ehf, and + energies from Eq. 
(p2|) for the two minima (oblate and prolate) of the two interactions, Sk3 and SG2. 
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Table 3. Results obtained from bare 2qp, TDA, and RPA calculations corresponding to the oblate and prolate 
solutions of SG2 and to the prolate solution of Sk3 in 74 Kr. The third (fourth) column contains the GT strength 
in units of [g\/ATi\ of the /3~ (/3 + ) decay summed up to E cut — 30 MeV. The fifth column contains the Ikeda sum 
rule obtained for this energy cut and the last column contains the energy-weighted summed strength in units of 
[MeV g\ /Air] corresponding to the (3 + decay. In all cases the results are for K = and K = 1 (within parentheses) . 
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Table 4. Total half-lives T 1/2 [minutes] of 74 Kr 
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